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ABSTRACT 


Context. Gravitational waves from black-hole (BH) merging events have revealed a population of extra-galactic BHs residing in short-period 
binaries with masses that are higher than expected based on most stellar evolution models — and also higher than known stellar-origin black holes 
in our Galaxy. It has been proposed that those high-mass BHs are the remnants of massive metal-poor stars. 

Aims. Gaia astrometry is expected to uncover many Galactic wide-binary systems containing dormant BHs, which may not have been detected 
before. The study of this population will provide new information on the BH-mass distribution in binaries and shed light on their formation 
mechanisms and progenitors. 

Methods. As part of the validation efforts in preparation for the fourth Gaia data release (DR4), we analysed the preliminary astrometric binary 
solutions, obtained by the Gaia Non-Single Star pipeline, to verify their significance and to minimise false-detection rates in high-mass-function 
orbital solutions. 

Results. The astrometric binary solution of one source, Gaia BH3, implies the presence of a 32.70 + 0.82 Mo BH in a binary system with a period 
of 11.6 yr. Gaia radial velocities independently validate the astrometric orbit. Broad-band photometric and spectroscopic data show that the visible 
component is an old, very metal-poor giant of the Galactic halo, at a distance of 590 pc. 

Conclusions. The BH in the Gaia BH3 system is more massive than any other Galactic stellar-origin BH known thus far. The low metallicity 
of the star companion supports the scenario that metal-poor massive stars are progenitors of the high-mass BHs detected by gravitational-wave 
telescopes. The Galactic orbit of the system and its metallicity indicate that it might belong to the Sequoia halo substructure. Alternatively, and 
more plausibly, it could belong to the ED-2 stream, which likely originated from a globular cluster that had been disrupted by the Milky Way. 


Key words. astrometry — binaries: spectroscopic — stars: black holes — stars: evolution — stars: massive — stars: Population П 


1. Introduction Stellar evolution models have difficulties in explaining such 


Since the first event detected in 2015 (Abbott et al. 2016) by 
the LIGO/Virgo collaboration, the detection of black-hole (BH) 
mergers via gravitational waves has uncovered the existence of a 
population of BHs residing in short-period binaries with masses 
AP than 30 Mo, ranging up to 85 Мо (Abbott et al. 20206, 


* Full Table B.1 and Table B.2 with Gaia epoch data are available at 
the CDS via anonymous ftp to cdsarc.cds.unistra.fr (130.79. 
128.5) or via https://cdsarc.cds.unistra.fr/viz-bin/cat/ 
J/A+A/686/L2 

* NASA Hubble Fellow. 

* Deceased. 


L2, page 2 of 23 


large masses for BHs of stellar origin: stars with an initial mass 
larger than 30 Мо are predicted to lose most of their mass during 
their evolution, due to the onset of strong winds, producing BHs 
with masses below 20 Мо (Vink 2008; Belczynski et al. 2007; 
Sukhbold et al. 2016). 

The masses of the merging BHs detected via gravitational 
waves are also larger than any known stellar-origin BHs in our 
Galaxy: all confirmed or candidate BHs of stellar origin in the 
Milky Way have typical masses around or below 10 Мо, with 
Сур X-1 (^20 Mo, Miller-Jones et al. 2021) being the most mas- 
sive one known thus far. However, the known stellar-origin BHs, 
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Table 1. Basic properties of Gaia BH3 from the Gaia DR3 catalogue. 


Parameter Value 

a [deg] 294.8278625082 (? 
Ta [mas] 0.051 

ó [deg] 14.9309796086 ( 
oe [mas] 0.052 

w [mas] 1.679 + 0.069 ®© 
џи“, [mas угт !] —22.235 + 0.062 
Ms [mas yr! —155.276 + 0.059 


G [mag] 11.2311 + 0.0028 
Ger — Grp [mag] 1.2156 + 0.0048 
Gays [mag] 10.2289 + 0.0122 
Ter [K] 5540755 

log g 3.08753 
[M/H] 16 
[o/Fe] 0.54 + 0.06 
RV [kms~!] —3332 + 3.4 


Notes. Astrophysical parameters are from GSP-Spec ANN. “The ref- 
erence epoch of DR3 coordinates is J2016.0 (JD 2457389.0). ФА zero- 
point correction (Lindegren et al. 2021a) of 35.4 uas has been applied 
to the parallax value given in the catalogue. 


mainly limited to short-period X-ray binaries, are only a very 
tiny fraction of the expected number of BHs in our Galaxy (~10° 
e.g. Olejak et al. 2020). In fact, stellar-origin BHs are hard to 
detect because most of them do not interact with a compan- 
ion. So, the lack of data on BHs with masses larger than 20 Mo 
could be due to an observational bias. Indirect methods such 
as gravitational microlensing have also yielded only one robust 
discovery of a single BH with mass, of about 7 Мо (Lam et al. 
2022; Sahu et al. 2022). 

The detection of mergers of BHs with masses larger than 
30 Мо can be reconciled with stellar evolution models if the 
progenitors of the high-mass BHs are low-metallicity stars 
(Mapelli et al. 2009; Belczynski et al. 2010a, 2016; Ziosi et al. 
2014; Fryeretal. 2012). The lack of metals substantially 
decreases the mass loss during the stellar lifetime (Vink 2008) 
and reduces the radius of the evolving progenitors (Hurley et al. 
2000; Belczynski et al. 2010a), the latter effect decreasing the 
probability of merging during the common-envelope phase 
(Belczynski et al. 2007) in binary systems. Finally, the higher 
mass of the BHs produced by low-metallicity progenitors is 
expected to decrease substantially or eliminate the natal kick 
strength at the birth of the BHs, preserving the binary as a bound 
system (Belczynski et al. 2010b). The maximum metallicity for 
the formation of the high-mass BHs is a matter of active debate, 
with some models predicting the formation of 30 Мо BHs even 
at solar metallicities (Bavera et al. 2023). 

Since its Data Release 3 (DR3, Gaia Collaboration 2023b), 
the Gaia mission (Gaia Collaboration 2016) has increased the 
number of detected stellar binary systems by two orders of 
magnitude (Gaia Collaboration 2023a; Halbwachs et al. 2023; 
Gosset et al. 2024) This has opened up the possibility of detect- 
ing BHs in binary systems that do not interact with their com- 
panion (see also Giesers et al. 2018; Saracino et al. 2022). More- 
over, the ability of Gaia to measure the astrometric orbit of such 
systems allows the measurement of the inclination of the orbit, 
providing a robust estimate of the mass of the dark companion. 
Two dormant BHs, Gaia BH1 and BH2 (El-Badry et al. 2023a,b; 
Tanikawa et al. 2023; Chakrabarti et al. 2023) were discovered 
in Gaia binaries of DR3. Gaia Data Release 4 (DR4) is expected 


to contain a larger number of binary systems than Gaia DR3; 
consequently, this will provide a greater number of BH-hosting 
systems, which will help to shed light on the BH population and 
the mechanisms in action in the BHs' formation. 

In this Letter, we report the serendipitous discovery of a 
nearby (~590 pc) binary system composed of an old, very metal- 
poor!, giant star orbiting a BH in 11.6 yr. The estimated BH 
mass, 33 Mo, is substantially higher than all known Galactic BHs 
and is in the mass range of the extra-galactic BHs detected by 
gravitational waves. 

The system was identified while validating the preliminary 
Gaia astrometric binary solutions produced in preparation for 
DR4 and subsequently confirmed by Gaia RVS radial-velocity 
data. We took the exceptional step of the publication of this paper 
based on preliminary data ahead of the official DR4 due to the 
unique nature of the discovery, which we believe should not be 
kept from the scientific community until the next release. An 
early disclosure will also enable an early and extensive follow- 
up by the community. 


2. Observations and analysis 
2.1. Properties of the source 


Gaia DR3 4318465066420528000 (also known as LS II +14 13 
and 2MASS J19391872+1455542), hereafter denoted as Gaia 
BH3, is a bright source in the constellation Aquila, known 
to be a high proper-motion star (Lépine & Shara 2005). Its 
basic properties from the Gaia DR3 archive are reported in 
Table 1. Its absolute magnitude and color (Riello et al. 2021; 
Sartoretti et al. 2023) identify it as a star climbing the giant 
branch (see Fig. 1). The source was analysed by the Astrophys- 
ical Parameters Inference System (Apsis, Creevey et al. 2023; 
Fouesneau et al. 2023). It has been identified as a G spectral- 
type star by the ESP-ELS algorithm (Sect. 11.3.7 of the online 
documentation, Ulla et al. 2022) and the GSP-Spec ANN param- 
eters (Recio-Blanco etal. 2023) indicate it as a metal-poor 
giant. No GSP-Phot result (Andrae et al. 2023) is published 
in the Gaia archive, while the parameters provided by GSP- 
Spec MatisseGauguin (Recio-Blanco et al. 2023) carry large 
uncertainties. 

The source is not known as a variable star in the litera- 
ture, neither in the AAVSO International Database, nor in the 
ASAS-SN database. We inspected ASAS-SN, ZTF, and TESS 
photometry, finding that the source does not present any sig- 
nificant periodic variability. The source was not observed with 
XMM-Newton, Chandra nor GALEX, nor it is present in the 
RAVE, APOGEE, LAMOST, or GALAH spectroscopic surveys. 
No eROSITA data have been made available yet for Gaia BH3, 
which belongs to the eastern Galactic hemisphere. 


2.2. Astrometry and orbital solution 


The system was identified while validating astrometric bina- 
ries orbital solutions produced by the Non-Single Star (NSS) 
pipeline in a preliminary run (identified as NSS 4.1), done in 
preparation of Gaia DR4. The NSS pipeline used in the NSS 4.1 
run is similar to the one used in the Gaia DR3, which is described 
in Halbwachs et al. (2023) and in Sect. 7.2.2 of the DR3 NSS 
documentation (Pourbaix et al. 2022), with improvements in the 
filtering of spurious solutions. The NSS pipeline processed astro- 
metric data produced by preliminary runs of the Intermediate 


! We use the nomenclature from Beers & Christlieb (2005) where very 
metal-poor are those stars having [Fe/H] < —2. 
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| е Gaia BH1 
a Gaia BH2 
Ф Gaia ВНЗ 


0.5 1.0 1.5 


(Ger - Grp)o 


Fig. 1. Gaia BH3 position in the Gaia color-magnitude diagram, com- 
pared with the position of Gaia BH1, BH2 and the low extinction 
(Ао < 0.05 mag) Gaia ОКЗ color-magnitude diagram. All extinctions 
are estimated through the Lallement et al. (2022) extinction map. 


2.0 2.5 


Data Update (IDU; see Fabricius et al. 2016) and the Astromet- 
ric Global Iterative Solution (AGIS; see Lindegren et al. 2021b) 
pipelines, covering the time range from JD 2456941.6218 to JD 
2458869.4177 (TCB?), for a total of about 64 months. The NSS 
4.1 run was executed on an input list of 10.4 million sources, 
chosen to be brighter than Grys = 14mag, and produced almost 
1.5 million orbital solutions. We note that the final NSS run for 
DR4 will extend to fainter magnitudes, and it is expected to pro- 
duce a much larger number of binary solutions. Further details 
on the NSS 4.1 run can be found in Appendix A. 

For each orbital solution, we computed the astrometric mass 
function from the angular semi-major axis of the photocentre 
orbit (ao), period (P), and parallax (w) as: 


fu = (ao/@)°(1 yr/ P Мо. (1) 


For an astrometric binary, the mass function depends on 
the masses of the components (Mi, М) and on their flux 
ratio, 72/71 (Halbwachs et al. 2023). For an invisible compan- 
ion (95/71 = 0) the mass function simplifies to: 
М» J 


———— p 
Mı + М Q) 


fu = М | 
from which it follows that М» > fw. Given fm and an estimate 
of the mass of the visible component (Mj), Eq. (2) can be used 
to solve for the mass of the dark companion (M2). 

Among the 1.5 million orbital solutions, Gaia BH3 yielded 
the largest mass function, 32.03 + 0.64 Mo, with a significance 
(а0/са,) of 48.1; no other solution has a mass function larger 
than 20 Me. In Fig. 2, we show the orbital solution and the resid- 
uals, from which the strength of the astrometric signal of the 
orbit, along with the robustness and quality of the solution can 


? TCB: Barycentric Coordinate Time, the time scale used here for all 
Gaia dates. 
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be appreciated. The Campbell orbital elements of the source are 
reported in the central column of Table 2. We note that the NSS 
pipeline used in this preliminary run produces Thiele-Innes ele- 
ments; the Campbell elements and their uncertainties were com- 
puted using the equations in Appendix A of Halbwachs et al. 
(2023). The astrometric mass function value and its uncertainty 
were computed using Monte Carlo resampling of the Thiele- 
Innes elements, the parallax and the period, in order to take 
into account the correlations between parameters; in particular, 
between ао and the period. To make sure this procedure would 
yield reliable results, we first checked that the correlations are 
sufficiently well behaved to allow for Monte Carlo resampling, 
following Sect. 6.1 of Babusiaux et al. (2023). 

A word of caution is necessary on the parallax value, and 
thus on the astrometric mass function: as in previous releases, the 
Gaia parallaxes are affected by a small bias (see Lindegren et al. 
2021a), but we do not have enough information at this stage to 
quantify the bias for the preliminary NSS solutions. As a conse- 
quence, the uncertainty on the mass function reported in Table 2 
is underestimated. 


2.3. Spectroscopy and combined orbital solution 


The Gaia RVS (Cropper et al. 2018) data of sources with an 
orbital solution from the NSS 4.1 run were processed with the 
DR4 operational RVS pipeline. Improvements with respect to the 
DR3 version will be described in the Gaia DRA documentation; 
the DR3 RVS pipeline is described in Sartoretti et al. (2018), 
Katz et al. (2023), and Sartoretti et al. (2022). It is worth men- 
tioning that the DR4 RVS pipeline includes the correction for the 
effect of the astrometric orbital motion, discussed in Holl et al. 
(2023); in particular, Sect. 3.3.1. The DR4 RVS data cover the 
time range from JD 2456863.9385 to JD 2458869.4177, namely, 
about 67 months. The pipeline produced 17 valid epoch radial 
velocities for Gaia BH3, reported in Appendix B, using a tem- 
plate spectrum with Т.ж = 6000 K, logg = 3.5 and [Fe/H] = 
—1.5 (see Blomme et al. 2017). The template parameters were 
estimated by the pipeline from the RVS spectrum itself. 

We used the DR3 NSS pipeline code to compute an orbital 
solution combining astrometric data and Gaia RVS radial 
velocities (nss. solution type = AstroSpectroSB1). The 
details of the adopted model are described in Sect. 7.7.3 of 
Pourbaix et al. (2022). We recall that in the combined solution 
model only the period, eccentricity, and periastron time are in 
common between the astrometric and the spectroscopic part of 
the solution. There is no constraint to impose that the semi-major 
axis of the photocentre orbit in AU (= ao/«) would be equal to 
the semi-major axis of the spectroscopic orbit o, as expected in 
the case of a dark companion. The consistency between ао and 
a, allows us to check whether the flux ratio is indeed compatible 
with a value of zero. 

As discussed in Sect. 2.2, the parallax derived by the NSS 
pipeline for the combined solution (w = 1.6808 + 0.0086 mas) 
is affected by a bias which we cannot quantify. In order to avoid 
its effect on the mass function, we estimate the latter using а 
instead of ао/о, namely: 


А ) um (3) 
1AU/ \ P SZ 


resulting in a value of 31.23 + 0.81 Мо. Assuming the equality 
between the photocentre and the spectroscopic orbit, we can also 
provide an alternative estimation of the parallax as 


| 


@=do/a, (4) 
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Fig. 2. Astrometric data of Gaia BH3. Top-left panel: Motion on the sky of the photocentre of the source, as seen by Gaia in the different CCD 
transits (dots), compared with the best fitting single-star solution from AGIS and the astrometric-binary solution from the NSS pipeline; the arrow 
indicates the direction of the proper motion. Bottom-left panel: Derived astrometric orbit of the photocentre, after a subtraction of parallax and 
proper motion, compared with the astrometric measurements. We note that only one-dimensional (1D) along-scan (AL) astrometry was used by the 
NSS pipeline. The position of the photocentre on the sky corresponding to each measurement is derived combining the measured one-dimensional 
AL position and the assumed orbital solution. The + signs show the barycentre and the position of the periastron, the dotted line shows the line of 
nodes, and the arrow indicates the direction of the motion along the orbit. In the top-right and bottom-right panels, we can see the residuals of the 
along-scan (AL) astrometric measurements for, respectively, the single-star solution and the binary-star solution. The vertical dot-dashed line in 


the bottom-right panel marks the time of the periastron passage. 


which results in a value of 1.6933 + 0.0164 mas. 

The Campbell orbital elements of the combined solution for 
Gaia BH3 are reported in Table 2. The combined solution is very 
similar to the astrometric solution, with slightly smaller uncer- 
tainties, and a better goodness-of-fit (GoF, in the HIPPARCOS 
sense, see Pourbaix et al. 2022), as a result of a stronger filter- 
ing of outliers. Radial velocities predicted by the combined solu- 
tion are compared with the measurements in Fig. 3. In Fig. 4, we 
show the combined and normalised RVS spectrum (see Seabroke 
et al. in prep.) compared with the template spectrum. 

Given the extreme value of the mass function of the sys- 
tem and the importance of its detection, a confirmation with 
ground-based observations was indispensable to discard the pos- 
sibility of a spurious solution. We thus observed Gaia BH3 
with the HERMES spectrograph (Raskin et al. 2011) mounted 


on the 1.2-m Mercator telescope at the Roque de los Mucha- 
chos Observatory (Spain), and with the SOPHIE spectrograph 
(Perruchot et al. 2008) mounted on the 1.93-m telescope of the 
Observatoire de Haute-Provence (France). A search in the ESO 
archive revealed that the source was observed with the UVES 
spectrograph (Dekker et al. 2000) mounted on the VLT. Details 
on the data reduction of these observations can be found in 
Appendix C. The spectra do not show any sign of the presence 
of a second component, nor of continuum filling of absorption 
lines; furthermore, no emission line was detected. Radial veloc- 
ities were derived for each ground-based observation and their 
values are reported in Table C.1. Although these values were not 
used to derive the orbital solution described above, they are in 
agreement with the predicted radial velocity within 0.5 km s^! ; 
this is less than the uncertainty of the orbital solution, as can be 
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Table 2. Campbell orbital elements of the Gaia BH3 system and the 
astrometric parameters of its barycentre. 


Parameter Astrometric Combined 
solution solution 

a [deg] 294.8278502411  294.8278502301 © 

Ta [mas] 0.060 0.054 

6 [deg] 14.9309190720 (? 14.9309190869 © 

с [mas] 0.086 0.074 

w [mas] 1.6747 + 0.0094 1.6933 + 0.0164 ®© 

и [mas yr !] —28.372 + 0.077 —28.317 x 0.067 

Ms [mas yr! —155.150 + 0.129 —155.221 € 0.111 

P [days] 4194.7 + 112.3 4253.1 + 98.5 

e 0.7262 + 0.0056 0.7291 + 0.0048 

ao [mas] 27.07 + 0.56 27.39 + 0.49 

i [deg] 110.659 + 0.107 110.580 + 0.095 

Т, (JD, TCB] 2458177.28 + 0.08 2458177.39 + 0.88 

Q [deg] 136.200 + 0.147 136.236 + 0.128 

w [deg] 7777 + 0.66 77.34 + 0.76 

a, [AU] Reg 16.17 + 0.27 

y[kms7!] ep —357.31 + 0.44 

Jm [Mo] 32.03 + 0.64 31.23 + 0.81 

GoF 2.17 —0.53 


Notes. (Тһе reference epoch of DR4 coordinates is J2017.5 (JD 
2457936.875). Derived as ao/a;, see text. 
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Fig. 3. Radial-velocity evolution of Gaia BH3. Top panel: Compari- 
son between the radial-velocity evolution predicted from the combined 
Gaia astrometric-spectroscopic binary model (blue solid line) and the 
epoch radial velocities measured with the Gaia RVS instrument (black 
filled circles), and ground-based measurements for Gaia BH3. Bottom 
panel: Radial-velocity residuals with respect to the binary solution com- 
pared with the 1-c- uncertainty of the predicted radial-velocity evolution 
(blue shaded area). The vertical dot-dashed line in both panels marks the 
time of the periastron passage. 


seen in Fig. 3. This result confirms the reality and accuracy of 
the orbital solution derived from Gaia data. 


2.4. Stellar parameters, abundances and Galactic orbit 


We derived new stellar parameters of the luminous compo- 
nent, using the Gaia DR3 photometry (G magnitude and 
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Fig. 4. Gaia RVS combined spectrum of Gaia BH3, in restframe, com- 
pared with the template spectrum. 


Table 3. Stellar parameters of Gaia BH3 derived in this work. 


Parameter Value 

Ter [K] 5212 + 80 
logg 2.929 + 0.003 
[Fe/H] —2.56 + 0.11 
[a/Fe] 0.43 + 0.12 
[M/H] —2.21 € 0.15 
log(L, / Lo) 1.208 + 0.030 
M, [Mo] 0.76 + 0.05 
R, [Rol 4.936 + 0.016 
Мсо [mag] 1.778 + 0.082 
(GBP = Grp)o [mag] 0.921 + 0.031 


Gpp — Grp colour), the parallax from the combined astrometric- 
spectroscopic solution, and the extinction (Ag) derived from 
the dust extinction maps of Vergely et al. (2022), with an iter- 
ative procedure described in Appendix D. The UVES spec- 
trum was used to determine the metallicity and abundances (see 
Appendix E). We then compared the extinction-corrected abso- 
lute G magnitude (Mao) and the dereddened colour (Свр — 
Ggp)o with the ones given by the isochrones libraries PARSEC 
(Bressan et al. 2012) and BaSTI (Pietrinferni et al. 2021). Thus, 
we derived the mass (M) of the visible component. The param- 
eters are reported in Table 3. 

In Fig. 5, we compare the expected spectral energy distri- 
bution (SED) with the Gaia XP spectrum (Carrasco et al. 2021; 
De Angeli et al. 2023; Montegriffo et al. 2023) and 2MASS pho- 
tometry. The agreement between the predicted SED and the Gaia 
XP spectrum is very good, with the only exception of the blue 
edge, where the XP spectrum is noisier. 

The abundances of Gaia BH3 (reported in Appendix E) show 
that the star is a-enhanced, as expected for a very metal-poor 
star. There is no trace of PC in the spectrum and the [Ba/Fe] 
is nearly solar, indicating the star has not been enriched by 
material processed in the CNO cycle, as expected if it had, for 
instance, accreted material from a companion star in the AGB 
phase. The star has no chemical peculiarity, except an enhance- 
ment of Eu ([Eu/Fe] = 0.52). Thus, it can be classified as 
an r-I neutron-capture-rich star, following the classification of 
Beers & Christlieb (2005). 
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Fig. 5. Gaia BH3 modelled SED, compared with the Gaia XP spec- 
trum and 2MASS photometry. The thin black line shows the unreddened 
model, while the thick line shows the SED assuming A, = 0.71 mag. 


Using the systemic radial velocity, proper motion, position, 
and distance, and assuming the Milky Way gravitational poten- 
tial from McMillan (2017), we may find that the source has 
a high-energy retrograde orbit? (Е = —1.29 x 10° km? 572, 
L, = -2.3 x 10? Крс kms^!, Lı = 1.05 x 10? kpc km 57!) 
in the Galaxy. Its kinematic characteristics are consistent with 
those of the halo substructure known as Sequoia (Myeong et al. 
2019), but are also in agreement to those of the recently dis- 
covered ED-2 stream (Dodd et al. 2023; Balbinot et al. 2023), 
a likely remnant of a globular cluster that was disrupted by the 
Milky Way. The metallicity of Gaia BH3 is more consistent with 
ED-2 ([Fe/H] = —2.6 + 0.2) than with the median metallicity of 
Sequoia ([Fe/H] ~ —1.7). 


3. Discussion 


With an estimated mass of 0.76 +0.05 Мо for the luminous com- 
panion, we derived a mass of: 


Mpy = 32.70 + 0.82 Mo , (5) 


for the dark companion. The observed luminosity of Gaia BH3 
is too low by several orders of magnitude to be compatible with 
the hypothesis that the companion is a main sequence star or 
even two main sequence stars in a close orbit. The estimated 
mass is also too large for one neutron star or two neutron stars 
in a close orbit, so we are left with the possibility of: (1) a single 
BH; (2) an inner binary containing two BHs; or (3) a BH and 
another compact object. Although the single BH is the simplest 
explanation, the hypothesis of an inner binary of two BHs can- 
not be excluded. Hayashi et al. (2023) proposed a method to test 
this hypothesis by detecting radial-velocity perturbations at the 
periastron. Using the Hayashi et al. (2023) formulation, we esti- 
mated radial-velocity perturbations* with a maximum amplitude 
of the order of 0.2 km s~!. Such perturbations are too small to be 
detected in the Gaia RVS data, but can be verified with ground- 
based instruments (see Nagarajan et al. 2024, for an application 


? Here we use the same conventions as in Myeong et al. (2019). 

^ We use Kaell — e) 7? as perturbations level estimation, (see Sect. 
2.1 in Hayashi et al. 2023), assuming an inner equal-mass binary with 
a circular coplanar orbit, and a period corresponding to the maxi- 
mum allowed by dynamic stability (126 days, according to Eq. (6) in 
Hayashi et al. 2023). 


to Gaia BH1). For the purposes of the subsequent discussion, we 
have adopted the single BH hypothesis as the most likely expla- 
nation. 

The estimated mass of the BH in Gaia BH3 makes it the most 
massive BH of stellar origin discovered in our Galaxy. It is strik- 
ing that the only BH with a mass larger than 20 Mọ found in the 
Gaia data so far is in orbit with a very metal-poor star, while such 
stars make up only a tiny fraction of the stars analysed in the NSS 
pipeline run (0.496 of sources which produced a binary solution 
have [М/Н] « —2 from DR3 GSP-Phot). Such stars also make 
up a small fraction of our Galactic halo (less than 5% accord- 
ing to Bonifacio et al. 2021) where this star and the majority of 
metal-poor stars are located. Although we can not exclude that 
this BH is the result of the merger of two less massive BHs, this 
discovery strongly supports the scenario where high-mass BHs 
are remnants of low-metallicity stars. The above considerations 
also raise the question of the maximum metallicity value for the 
formation of high-mass BHs, which in Belczynski et al. (2016) 
is identified at [M/H] = —1. The much lower metallicity of Gaia 
BH3 may be an indication that high-mass BHs form only at very 
low metallicities rather than at moderately low ones. 

An in-depth discussion of the possible formation scenarios 
for this binary system is beyond the scope of the paper; never- 
theless, a few aspects ought to be highlighted. As discussed in 
El-Badry et al. (2023b,a), the formation of the Gaia ВНІ and 
BH2 systems as isolated binaries is unlikely. This is also true for 
the recently discovered Gaia NS1 system (El-Badry et al. 2024), 
composed of a high-mass neutron star and a low-metallicity 
star. Given the size of their orbits, these systems should have 
experienced a common-envelope phase and then a mass trans- 
fer toward the light companion, which would then have resulted 
in much closer orbits than the observed ones. For Gaia BH2, the 
common-envelope phase could have been avoided if the BH pro- 
genitor was more massive than 65 Mo. In the case of Gaia BH3, 
the present-day minimum separation is of the order of 1000 Ro 
and the common-envelope phase could not have been avoided 
because models predict that the BH progenitor becomes a red 
supergiant even at 150 Mọ (Chen et al. 2015). Similarly to Gaia 
ВНІ and BH2, the chemical composition of the luminous com- 
ponent does not show any unusual abundance; in particular, the 
absence of ІС and the observed [Ba/Fe] point toward a lack of 
contamination by the BH progenitor during its evolution. The 
observed enhanced Eu abundance could be due to the contam- 
ination from the SN at the birth of the BH, but also due to the 
medium in which the star formed. An alternative formation sce- 
nario, proposed to explain the Gaia BH1 and BH2 systems, is 
that the BH acquired the low-mass companion via dynamical 
exchange in a dense environment (see for example Rastello et al. 
2023; Tanikawa et al. 2024; Di Carlo et al. 2024). Such a sce- 
nario might be supported by the probable association of Gaia 
BH3 with the ED-2 stream, which could be a remnant of a glob- 
ular cluster (Dodd et al. 2023; Balbinot et al. 2023). 


4. Conclusions 


In this Letter, we present the discovery of a wide binary com- 
posed of a very metal-poor giant orbiting a dark object of 33 Mo, 
using Gaia preliminary DR4 astrometric data, corroborated by 
Gaia spectroscopy. Most probably, the massive dark object is 
a single black hole (BH). The 33 Mọ of the BH mass makes it 
the most massive BH of stellar origin discovered in our Galaxy. 
All Galactic BHs that reside in short-period X-ray binaries have 
masses generally below 10 Мо (e.g. Corral-Santana et al. 2016), 
except Cyg-X1 (Мвн ~ 20 Mo). Even the first two dormant BHs 


L2, page 7 of 23 


Gaia Collaboration: A&A, 686, L2 (2024) 


discovered by Gaia in wide astrometric orbits have masses of 
about 10 Mo. The mass of Gaia BH3 puts it in the mass range 
of the BHs discovered by gravitational waves (e.g. Abbott et al. 
2023), and, in fact, it is close to the peak of the observed mass 
distribution for the merging BHs (e.g. Farah et al. 2024). The 
metallicity of the system supports the scenario (Belczynski et al. 
2016) that the high-mass BHs observed by LIGO/Virgo/KAGRA 
(Abbott et al. 20202) are the remnants of metal-poor stars. 

The discovered system, with its extremely low-mass ratio, 
wide orbit, and specific chemical composition, can also provide 
constraints for stellar evolution and binary models. As in the case 
of the Gaia BH1 and BH2 systems, the formation scenario as an 
isolated binary appears unlikely and alternative scenarios should 
be considered. The BHs discovered by Gaia in wide binaries in 
our Galaxy and those detected by LIGO/Virgo/K AGRA in exter- 
nal galaxies (i.e. BH merger events of extremely short-period 
binaries) constitute two ends of the BH population. When stud- 
ied together, they can help to formulate a comprehensive view of 
BH formation and the evolution of their progenitors. 

Finally, the bright magnitude of the system and its relatively 
small distance makes it an easy target for further observations 
and detailed analyses by the astronomical community. This dis- 
covery should be also seen as a preliminary teaser for the con- 
tent of Gaia DR4, which will undoubtedly reveal other binary 
systems hosting a BH. 
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Appendix A: Astrometric processing in the NSS 4.1 
run 


Here, we provide the details of the Non-Single Star (NSS) 
pipeline used in the preliminary run NSS 4.1, underlining this 
is not the final version of the NSS pipeline for the generation of 
Gaia DR4 data. 

The pipeline is similar to the one used in Gaia DR3 
(Halbwachs et al. 2023), with a few updates — the main dif- 
ferences are in the improvement of the filters to remove spu- 
rious solutions. In particular, two new filters were introduced 
to remove sources which are partially resolved: the first one 
filters out sources with significant scan-angle-dependent sig- 
nals (Holl et al. 2023) in the G flux before attempting an astro- 
metric solution; the second one removes solutions which have 
periods that correspond to combinations of the Gaia spacecraft 
precession frequency and the yearly frequency (see Holl et al. 
2023, Sect. 4.1). These new filters allow us to relax the fil- 
ters based on the significance and on @/o,, described in 
Halbwachs et al. (2023), which are replaced by the following 
criteria: the goodness-of-fit (GoF) must be smaller than 15, the 
eccentricity error smaller than 0.2, semi-major axis significance 
a0/Ta > 5, and 


c/c, > max (15, 208.02 - log(P/1 day) + 548.03) . (А.1) 


The NSS 4.1 run was executed using astrometric data pro- 
duced with preliminary runs (namely IDU 4.1 and AGIS 4.1) 
of the Intermediate Data Update (IDU; see Fabricius et al. 
2016) and the Astrometric Global Iterative Solution (AGIS; see 
Lindegren et al. 2021b) pipelines. The preliminary astrometry 
provided by AGIS 4.1 covers the entire range of DR4, i.e. from 
JD 2456863.9385 to JD 2458869.4177, for a total of about 67 
months, but only data after JD 2456941.6218 were used in the 
NSS 4.1 run. The local perspective effect (Halbwachs 2009) was 
not included in the model for the run NSS 4.1 and the variability- 
induced mover (VIM) solutions were not attempted. 

The NSS 4.1 run was executed on a list of 10 450 939 sources 
chosen with the following criteria: the source must be brighter 
than Grys = 14mag and either G < 18mag, an astromet- 
ric renormalised unit weight error (RUWE) larger than 1.05, or 
w > 5mas and RUWE > 0.9, and a number of visibility peri- 
ods used in AGIS solution larger than 11. In order to exclude 
partially resolved sources, sources with a percent of successful 
Image Parameter Determination (IPD) windows with more than 
one peak larger than 10% or with amplitudes of the IPD GoF ver- 
sus the scan angle larger than 0.2, were excluded (see Gaia DR3 
archive documentation for details on the above quantities). The 
run produced a total of 1 469 196 orbital solutions. 

The selection function of the NSS 4.1 run is not trivial to 
characterise, thus, it is not possible to estimate how common or 
rare are systems like Gaia ВНЗ. If we push Gaia ВНЗ to the dis- 
tance corresponding to the cut in magnitude (Grys = 14 mag 
would correspond to a distance of 3.3 kpc, ignoring the extinc- 
tion), the semi-major axis of the orbit would be 4.8 mas, which is 
still very large with respect to the precision of Gaia epoch mea- 
surements at that magnitude. However, given that the DR4 time 
range covers only half of the orbital period, the resulting signif- 
icance could be below the acceptance thresholds. We note that 
during the Gaia DR3 preparation, Gaia BH3 produced an astro- 
metric acceleration solution (Halbwachs et al. 2023) and an SB1 
solution (Gosset et al. 2024), which were both discarded from 
the release due to a low significance, because the orbital period 
was much longer than the Gaia DR3 time span and the periastron 
passage was not covered by the DR3 time range. 
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If we consider a Gaia BH3-like system but with an orbital 
period similar or shorter of the DR4 time range (i.e. below 2000 
days), the significance would be always higher than the accep- 
tance threshold solution, with the exception of very short («20 
day) periods. We note that the NSS 4.1 is limited to periods 
larger than 10 days. For binaries with shorter periods, the giant 
would almost fill its Roche lobe and the source would probably 
be detected as a X-ray source. 

Although it has not yet been finalised, the input list for 
DR4 will be significantly larger than for NSS 4.1, probably 
built as the sum of a volume-limited sample and a G < 18 
magnitude-limited sample, as indicated above, though without 
the Grys < 14 mag criterion. The motivation for the dedicated 
NSS 4.1 run and the reason for the latter criterion was the analy- 
sis of the effect of a deviation of the astrometry from the assumed 
single-star model on the calibration of the spectroscopic instru- 
ment. 


Appendix B: Gaia epoch data 


Here, we describe the Gaia epoch astrometric data and epoch 
radial velocities used to produce the binary solution of Gaia 
BH3. 

The astrometric measurements of Gaia BH3 are provided in 
Table B.1. They were produced from preliminary pipelines, pro- 
visional instrument models and calibrations; as a consequence 
they will not be identical (but still similar) to the correspond- 
ing data to be produced and published for this star with DR4. 
Furthermore, the final epoch astrometry table in DR4 will con- 
tain many additional details and quality diagnostics on the indi- 
vidual measurements. A full explanation of the epoch astrome- 
try is beyond the scope of the present short appendix. We refer 
the reader to the Gaia Technical Document Lindegren & Bastian 
(2022) and to Lindegren et al. (2012). 

Each Gaia epoch astrometry record in Table B.1 corresponds 
to a transit of the source on one of the CCDs of the AF instru- 
ment. The Table is arranged as follows. Col. 1: transit id?,a 
unique identifier assigned to each detected celestial light source 
as its image transits the Gaia focal plane; Col. 2: AF CCD strip; 
Col. 3: Barycentric time, in JD, corresponding to the middle of 
the 4.41-second CCD exposure time; Col. 4: Along-scan posi- 
tion of the photocentre, with its associated uncertainty, which 
corresponds to the longitude of the observed photocentre in a 
2D tangential coordinate system having its origin at a reference 
equatorial position, and having the axis of its longitude coordi- 
nate oriented corresponding to the scanning direction; Col. 5: 
Parallax factor, namely, the quantity by which it is necessary 
to multiply the parallax in order to obtain the contribution to 
the along-scan position due to the orbit of the spacecraft with 
respect to the Solar System barycentre; Col. 6: Scan angle, which 
is the angle of the scanning direction with respect to local ICRS 
North; Col. 7: Outlier flag, indicating whether the measurement 
was considered as an outlier (flag = 1) by the NSS pipeline and 
filtered out, or not (flag = 0), when solving for the astrometric- 
spectroscopic combined solution. 

The reference position (ao, бо), in 
Lindegren & Bastian (2022), for Gaia BH3 is 


the sense of 


ао = 294°82784900557243 , бо = 14°930918410309376 , (B.1) 


while the reference time is J2017.5 (JD 2457936.875). 


5 A decoder for the transit_id is available on-line at https:// 
gaia.esac.esa.int/decoder/transitidDecoder. jsp 
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transit_id AF Time Along-scan position Parallax factor Scan angle Outlier 
strip [BJD, TCB] [mas] [deg] flag 
20114916805338633 1 2456958.110978 147.066 + 0.370 0.70827985 —59.04672662 0 
20114916805338633 2 2456958.111034 146.696 + 0.231 0.7082799 1 —59.04677105 0 
20114916805338633 3 2456958.111091 146.685 + 0.183 0.70828003 —59.04681553 0 
20114916805338633 4 2456958.111156 146.557 + 0.151 0.70828015 —59.04686711 0 
20114916805338633 5 2456958.111203 146.396 x 0.097 0.70828021 —59.04690435 0 
20114916805338633 6 2456958.111259 146.374 + 0.088 0.70828027 -59.04694869 0 
20114916805338633 7 2456958.111325 146.436 + 0.100 0.70828038  —59.04700014 0 
20114916805338633 8 2456958.111372 146.060 + 0.100 0.70828044 -59.04703723 1 
20114916805338633 9 2456958.111437 146.256 + 0.130 0.70828056 —59.04708855 0 
20119009095238725 1 2456958.184993 146.244 + 0.364 0.70847327 . —59.11055165 0 
20119009095238725 2 2456958.185043 146.148 + 0.296 0.70847344 ^ —59.11059886 0 
20119009095238725 3 2456958.185105 146.192 + 0.267 0.70847368 —59.11065699 0 
20119009095238725 4 2456958.185162 146.271 + 0.225 0.70847386 | —59.11070963 0 
20119009095238725 5 2456958.185218 146.047 + 0.273 0.70847410 | —59.11076223 0 
20119009095238725 6 2456958.185274 146.068 + 0.158 0.70847428 —59.11081482 0 
20119009095238725 7 2456958.185330 146.023 + 0.177 0.70847452 —59.11086742 0 
20119009095238725 8 2456958.185386 0.70847470 0 


146.095 + 0.191 


—59.11092003 


Notes. The full table is available at the CDS. 


Table B.2. Gaia BH3 epoch radial velocities from Gaia RVS. 


transit id Time Radial velocity 
[BJD, TCB] [ кт 571] 

22989619581449144 2457010.096166 338.62 + 1.42 
22993711812243969 2457010.170169 341.08 + 1.70 
27565894819405940 2457092.856649 -338.73 + 1.41 
27569987068420645 2457092.930660 -340.63 + 1.43 
29624665178404300 2457130.091203  —33528 + 2.20 
37997278005744779 2457281.508798 | —330.88 + 1.24 
52687653750370648 2457547.177378 | —-320.11 + 1.61 
60072325242450363 2457680.722202  —319.09 + 1.87 
61702852954182073 2457710.207118 -315.11 + 1.12 
68202800816349412 2457827.754668 -305.82 + 1.53 
81052816639971536 2458060.142534  —290.29 + 1.42 
82486181452083612 2458086.062277 | —294.93 + 1.53 
82490273693622584 2458086.136278 | —-294.67 + 1.75 
100312802889225719 2458408.450948  —404.70 + 1.14 
111669627072959442  2458613.834346 -399.65 + 1.21 
112765863248582085 2458633.660613 —398.97 + 1.10 
121996985218154856 2458800.596176 -393.59 + 2.00 


Notes. The table is also available at the CDS. 


Epoch RVS radial velocities, reported in Table B.2, were 
produced with the final pipeline, but not finalised with the post- 
processing; their values or uncertainties may slightly differ from 
the final DR4 values. As for the astrometry, the final epoch radial 
velocity table in DR4 will contain additional details and quality 
diagnostics on the individual measurements. Each epoch radial 
velocity record in Table B.2 corresponds to a transit of the source 
on the RVS CCDs. The provided observation time of the radial 
velocity corresponds to the mean of the observation times of the 
three CCDs used to collect spectra in the RVS during the transit. 


A public code, illustrating the use of epoch astrometric and 
radial velocity data to produce an orbital solution for Gaia BH3 
is available online. 


Appendix C: Ground-based spectroscopy 


Table C.1. Gaia BH3 epoch radial velocities from ground-based obser- 
vations. 


Instrument Time Radial velocity 
[BJD, UTC] [km 571] 
UVES 2459158.5333  —383.20 + 0.30 
HERMES 24601434948  —364.05 + 0.20 
SOPHIE 2460192.4359 -363.15 + 0.02 
HERMES 2460193.4833  —36324 + 0.20 


We observed Gaia BH3 with the HERMES spectrograph 
(Raskin et al. 2011) mounted on the 1.2-meter Mercator tele- 
scope at the Roque de los Muchachos Observatory (Spain), at 
two dates (17 July 2023 and 7 September 2023), taking two con- 
secutive exposures of 2700s each night. The spectra have a spec- 
tral coverage from 377 to 900 nm, a resolving power of 85 000, 
and a S/N ~ 43 at 520 nm. 

We observed Gaia BH3 also with the SOPHIE spectrograph 
(Perruchot et al. 2008) mounted on the 1.93-meter telescope of 
the Observatoire de Haute-Provence (France) on 4 September 
2023. The source was observed with a single exposure of 6000 s; 
the spectra cover the range 387 to 694nm with a resolving 
power of 40000, and have a S/N ~ 66 at 520nm. HERMES 
and SOPHIE spectra of Gaia BH3 are available on request from 
the corresponding author. 

A search in the ESO archive revealed that the source was 
observed with the UVES spectrograph (Dekker et al. 2000) 


6 https://www.cosmos.esa. int/web/gaia/gaia-bhthree 
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Fig. C.1. Residual intensity of the Gaia BH3 UVES spectrum in the magnesium triplet region (top panel), in the Сап H+K region (middle panel), 


Hj (bottom-left) and Ha (bottom-right). 


mounted on the VLT, on 5 November 2020, in the program 
106.21JJ.001 proposed by T. Matsuno and collaborators. The 
aim of the program was to derive a complete chemical inven- 
tory of stars belonging to Galactic accretion events. The expo- 
sure time of the UVES spectrum was 900s in the 3904-580 set- 
ting (spectral coverage 326 to 454 nm and 476 to 684 nm) with 
a slit of 0.7", producing a resolving power of 58 000 in the blue 
and 62 000 in the red and a S/N — 100 at 520 nm. 

For HERMES and SOPHIE observations, radial velocities 
were derived by computing the cross-correlation functions with 
a G2 mask, while for the UVES spectrum the radial veloc- 
ity was derived via template matching. The barycentric radial- 
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velocity values are reported in Table C.1, and they are in good 
agreement with the orbit derived from Gaia data. The spec- 
tra do not show any sign of the presence of a second compo- 
nent, nor any emission line. The UVES normalised spectrum 
of the Gaia BH3 in selected spectral regions is shown in 
Fig. C.1. 


Appendix D: Derivation of stellar parameters 


Here we describe the iterative procedure used to derive the stellar 
parameters of the luminous component in Gaia BH3. The pro- 
cedure is similar to the one adopted in Lombardo et al. (2021). 
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Fig. D.1. Extinction in the direction of Gaia BH3, as function of the 
distance. The extinction was derived with the maps with correlation 
length of 10 pc (black solid line) and 25 pc (red dot-dashed line) from 
Vergely et al. (2022); the vertical shaded region shows the distance 
range of Gaia BH3. 


We computed the emerging flux for a grid in Teg, 1059 
and [Fe/H], of 1D plane-parallel model atmospheres, using 
ATLAS 9 (Castelli et al. 2003; Kurucz 2005). The fluxes were 
then converted to spectral energy distributions by multiply- 
ing them by a factor 47(Ro/10 pe)’, where Ro is the solar 
radius, as done in Casagrande & VandenBerg (2014). All mod- 
els were computed assuming [a@/Fe] = 0.4. For each model, we 
computed theoretical values of the color Свр – Сър, the bolomet- 
ric correction (BCg), and the extinction coefficients Ас /Ао and 
Е(Свр — Свр)/Ао, using the average Milky Way reddening law 
from Fitzpatrick et al. (2019). 

The iterative procedure starts with first-guess metallicity, 
temperature and gravity from Gaia DR3 values reported in 
Table 1, a first-guess mass OM, of 0.8 Mo, and a given redden- 
ing of Ао. We use the above parameters to obtain E(Gpp — Grp) 
from the grid, which is then used to obtain a dereddened Свр — 
Grp colour, (Свр — Ggp)o , the extinction Ас, and the bolomet- 
ric correction BCg. We then compare the (Gap — Сър)о value to 
theoretical colours in the grid to derive a new effective tempera- 
ture, and we use the Stefan-Boltzmann equation to derive a new 
surface gravity: 


(D.1) 


M, Tet 
1 = log —- +4log — + log go + 
ogg mu "lege + logge 


bo 
+ 0.4(G – Ae BCg lloc) +2 log 1000 mas ` 

The procedure is repeated to convergence in Тең and log g which 
is achieved after a few iterations. The parameters Te and log g 
are then used to derive the metallicity [Fe/H] as described in 
Sect. E and the process repeated. 

The Ter, logg, and [Fe/H] derived with the above proce- 
dure depend mainly on the choice of Ap. Gaia ВНЗ has a low 
Galactic latitude (b = —3.49?), located in a zone with a rel- 
atively high gradient of Ao, according to extinction maps of 
Vergely et al. (2022), as can be seen in Fig. D.1. Adopting a dis- 
tance of 590.6 + 5.8 pc, we obtain Ag = 0710100 mag for a 
correlation length of 10 pc and Ag = 0.666 x 0.017 mag for a cor- 
relation length of 25 pc. We thus adopted Ag = 0.71 + 0.07 mag. 

An updated value for M, can be then estimated by compar- 
ing the absolute G magnitude, Мв and the colour (Свр — Grp)o 
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Fig. D.2. Comparison between the position of Gaia BH3 and isochrones 
in the colour-magnitude diagram. The colours of the symbols (filled 
circles for BaSTI and filled squares for PARSEC), on the two isochrones 
sets (12 and 14 Gyr), correspond to the stellar masses. 


with the ones given by theoretical isochrones. The procedure for 
the determination of stellar parameters is finally repeated with 
the updated value of M,. 

The value of M, depends mainly on the assumed age 
and isochrone set, while it has a very small dependency 
on the assumed Ap. We used isochrones from both PAR- 
SEC’ (Bressan et al. 2012) and BaSTI? (Pietrinferni et al. 2021) 
libraries. For BaSTI, we adopt [Fe/H] = -2.5, [a/Fe] = 0.4 
and [M/H] = -2.18, while for PARSEC isochrones, which 
are only available with no a-enhancement, we use metallicity 
[Fe/H] = —2.18, i.e. we scale the [Fe/H] to match the [M/H] 
of BaSTI isochrones and take into account the contribution of 
a-elements to the total metallicity (e.g. Salaris et al. 1993). The 
comparison between isochrones and Gaia ВНЗ in the Mg ver- 
sus (Свр — Grp)o diagram is shown in Fig. D.2. 

The value of M, goes from 0.758 to 0.793 Mo for PARSEC 
and 0.723 to 0.755 Mo for BaSTI, for isochrone of ages 12 
and 14 Gyr, respectively. Younger ages would result in higher 
masses, but also bluer colours. We then estimate a mass M, = 
0.76 + 0.05 Mo as the mass for the visible companion. 

With the procedure described above, we obtain the follow- 
ing parameters: Tep = 5211 + 80K, logg = 2.929 + 0.003, 
[Fe/H] = —2.56 + 0.12. From the Мот and Сат abundances, 
we derived an a-enhancement of [a/Fe] = 0.43 + 0.12. Using 
the relation between iron content, enhancement and metallicity 
from Tantalo et al. (1998), we obtain [M/H] = -2.21 + 0.15. 

It can be seen that the source is slightly redder and cooler 
than what is predicted by the models, albeit not significantly. 
In order to check that the Te that we determined above is not 
underestimated, we derived an alternative temperature estima- 
tion from the excitation equilibrium of the Fer lines, including 
the NLTE corrections by Frebel et al. (2013). With this method, 
we obtained Te ~ 5100 K, which is even cooler, confirming that 
our estimation is not underestimated. A more detailed analysis of 
the stellar parameters is outside the scope of this work. 


7 http: //stev.oapd.inaf.it/cgi-bin/cmd 
8 http://basti-iac.oa-abruzzo.inaf.it/isocs.html 
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Table E.1. Abundances of Gaia BH3 from UVES spectrum. 


Ton Mines A(X) [X/H] TAX) NLTE 
(1) (2) corr. 

CI 6.11 -2.39 0.10 0.17 As 
Met 5 5.46 -2.08 0.05 0.04 0.067% 
Cal 20 413 -2.20 010 0.05 0.14? 
Sc Il 6 0.76 -2.34 0.11 0.03 

Til 17 2.67 -2.23 0.06 0.11 

Till 26 2.80 -2.10 0.08 0.04 

VI 3 1.34  -2.66 0.13 0.10 

Vil 6 161 -2.39 0.07 0.03 

Cri 5 2.86 -2.78 0.01 0.08 

Сгп 3 3.23 -2.41 0.06 0.03 

MnI 7 2.45 -2.92 0.04 0.04 

Mnn 1 2.79 -2.58 Se 0.08 ge 
Fel 223 496 -2.56 0.12 0.08 0.10° 
Fer 15 5.03 -2.449 0.3 0.01 к 
Col 7 2.60 -2.30 0.09 0.11 

Nil 18 3.68  -2.55 0.12 0.06 - 
ZnI 1 2.22 -2.40 a 0.02 0.124 
Sril 2 0.50 -2.42 0.05 0.08 8 
ҮП 10 -0.44 -2.65 0.16 0.05 

Zr Il 6 0.34 -2.28 0.12 0.05 а 
Вап 4 -0.28 -2.45 0.13 0.06 -0.10° 
Тап 4 -1.22 -2.36 0.07 0.05 

Eur 1 -1.45 -1.97 - 0.05 


Notes. The values A(X) are expressed in the form A(X) = log(X/H) + 
12, while [X/H] = log (А(Х)/А(Х),). The uncertainties on the abun- 
dances LO aux) are reported as (1) line-to-line dispersion or (2) effect of 
Teg uncertainty. NLTE corrections are from “Bergemann et al. (2017) 
’Mashonkina et al. (2007) ‘Bergemann et al. (2012) ^Sitnova et al. 
(2022) *Korotin et al. (2015) 


Appendix E: Abundances 


We estimated the metallicity and abundances of the luminous 
component of Gaia BH3 from the UVES spectrum, which is of 
higher quality than our HERMES and SOPHIE spectra. 

The abundances were derived with the code MyGIsFOS 
(Sbordone et al. 2014), with the exception of those for C, Ba 
and Eu; the carbon abundance was derived by fitting the Fraun- 
hofer G-band of the CH molecule, while line-profile fitting was 
used to determine Ba and Eu abundances. The resulting values 
are listed in Table E.1. In the table, we provide the line-to-line 
scatter and the variation of the abundance corresponding to the 
uncertainty on the effective temperature. We adopted the solar 
abundances for C, Fe and Eu from Caffau et al. (2011), and those 
from Lodders et al. (2009) for the other elements. We derived 
a value of 1.19 kms ^! for the micro-turbulence, by forcing the 
same Fe abundance from Fe I lines of different equivalent width. 

We sought the non-local thermal equilibrium (NLTE) correc- 
tion for Mg I and Car in Kovalev et al. (2018). For the four Mg I 
lines available, we derived a NLTE correction of 0.06. Eleven 
lines of Cat provided a NLTE correction of 0.15. The five Crit 
lines available in the database provide a large NLTE correction: 
0.46. For the Mn I features, the NLTE correction is 0.54. The 67 
Fe I lines available provide a mean NLTE correction of 0.1. The 
NLTE effect on the Znt line at 481 nm in Sitnova et al. (2022) 
and the Ba in Korotin et al. (2015) are small (see Table E.1). 

The [Fe/H] ratio we derived from the UVES spectrum 
([Fe/H] = —2.56 + 0.12 from 223 Fe I lines) 15 in perfect agree- 
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Fig. E.1. Residual intensity of the Gaia BH3 UVES spectrum (black 
line) compared with the modelled spectrum (thin red line) in the Eu 
region. 
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ment with the value derived from the SOPHIE spectrum ([Fe/H] 
= —2.57 +0.12 from 121 Fe 1 lines) and from the HERMES spec- 
trum ([Fe/H] = —2.54 + 0.14 from 156 Fe I lines). A good agree- 
ment was also obtained for the other elements, whose abundance 
is based on several lines. The star, as expected for the metal-poor 
regime, is enhanced in a elements. 

There is no trace of ?C in the spectrum, so the star has not 
been enriched by material processed in the CNO cycle, as it 
would if it had, for instance, accreted material from a companion 
star in the AGB phase. 

The star has no chemical peculiarity, except a slight enhance- 
ment in Eu: [Eu/Fe] = 0.52. When coupled with [Ba/Fe] = 
0.11, this classifies this star as an r-I neutron-capture-rich star, 
following the classification of Beers & Christlieb (2005). The 
UVES spectrum in the region of the Eu is shown in Fig. E.1. 

According to the stellar parameters, the star is expected 
to have a Li abundance on the Mucciarelli plateau (see 
Mucciarelli et al. 2012). At the wavelength of the Li feature at 
607 nm, we see an absorption line compatible with an abundance 
of A(Li) = 1.2, but the shape is not comparable with the synthetic 
spectrum, showing an absorption on the blue side of the feature. 
The Li feature is also visible in the SOPHIE spectrum, but with 
a lower S/N than in the UVES spectrum, while the S/N of the 
HERMES spectrum at 760 nm is too low. 
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Planck Institute for Astrophysics, Max Planck Institute for 
Extraterrestrial Physics, New Mexico State University, New 
York University, Ohio State University, Pennsylvania State 
University, University of Portsmouth, Princeton University, 
the Spanish Participation Group, University of Tokyo, Uni- 
versity of Utah, Vanderbilt University, University of Vir- 
ginia, University of Washington, and Yale University; 

— the thirteenth release of the Sloan Digital Sky Survey (SDSS 
DR13, Albareti et al. 2017). Funding for SDSS-IV has been 
provided by the Alfred P. Sloan Foundation, the United 
States Department of Energy Office of Science, and the Par- 
ticipating Institutions. SDSS-IV acknowledges support and 
resources from the Center for High-Performance Comput- 
ing at the University of Utah. The SDSS web site is https: 
/ ттт. sdss . org/. SDSS-IV is managed by the Astrophys- 
ical Research Consortium for the Participating Institutions of 
the SDSS Collaboration including the Brazilian Participation 
Group, the Carnegie Institution for Science, Carnegie Mel- 
lon University, the Chilean Participation Group, the French 
Participation Group, Harvard-Smithsonian Center for Astro- 
physics, Instituto de Astrofísica de Canarias, The Johns 
Hopkins University, Kavli Institute for the Physics and 
Mathematics of the Universe (IPMU) / University of 
Tokyo, the Korean Participation Group, Lawrence Berke- 
ley National Laboratory, Leibniz Institut für Astrophysik 
Potsdam (AIP), Max-Planck-Institut für Astronomie (MPIA 
Heidelberg), Max-Planck-Institut für Astrophysik (MPA 
Garching), Max-Planck-Institut für Extraterrestrische Physik 
(MPE), National Astronomical Observatories of China, New 
Mexico State University, New York University, University of 
Notre Dame, Observatário Nacional / MCTI, The Ohio State 
University, Pennsylvania State University, Shanghai Astro- 
nomical Observatory, United Kingdom Participation Group, 
Universidad Nacional Autónoma de México, University of 
Arizona, University of Colorado Boulder, University of 
Oxford, University of Portsmouth, University of Utah, Uni- 
versity of Virginia, University of Washington, University of 
Wisconsin, Vanderbilt University, and Yale University; 

— the second release of the SkyMapper catalogue (SkyMap- 
per DR2, Onkenetal. 2019, Digital Object Identifier 
10.25914/5ce60d31ce759). The national facility capability 
for SkyMapper has been funded through grant LE130100104 
from the Australian Research Council (ARC) Linkage 
Infrastructure, Equipment, and Facilities (LIEF) programme, 
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awarded to the University of Sydney, the Australian National 

University, Swinburne University of Technology, the Univer- 

sity of Queensland, the University of Western Australia, the 

University of Melbourne, Curtin University of Technology, 

Monash University, and the Australian Astronomical Obser- 

vatory. SkyMapper is owned and operated by The Australian 

National University's Research School of Astronomy and 

Astrophysics. The survey data were processed and provided 

by the SkyMapper Team at the Australian National Univer- 

sity. The SkyMapper node of the All-Sky Virtual Obser- 
vatory (ASVO) is hosted at the National Computational 

Infrastructure (NCI). Development and support the SkyMap- 

per node of the ASVO has been funded in part by Astron- 

omy Australia Limited (AAL) and the Australian Govern- 
ment through the Commonwealth's Education Investment 

Fund (EIF) and National Collaborative Research Infrastruc- 

ture Strategy (NCRIS), particularly the National eResearch 

Collaboration Tools and Resources (NeCTAR) and the Aus- 

tralian National Data Service Projects (ANDS); 

— Ше Gaia-ESO Public Spectroscopic Survey (GES, 
Gilmore et al. 2022; Randichetal. 2022). The Gaia- 
ESO Survey is based on data products from observations 
made with ESO Telescopes at the La Silla Paranal Observa- 
tory under programme ID 188.B-3002. Public data releases 
are available through the ESO Science Portal. The project 
has received funding from the Leverhulme Trust (project 
RPG-2012-541), the European Research Council (project 
ERC-2012-AdG 320360-Gaia-ESO-MW), and the Istituto 
Nazionale di Astrofisica, INAF (2012: CRA 1.05.01.09.16; 
2013: CRA 1.05.06.02.07). 

The GBOT programme (GBOT) uses observations collected 
at (1) the European Organisation for Astronomical Research in 
the Southern Hemisphere (ESO) with the VLT Survey Tele- 
scope (VST), under ESO programmes 092.B-0165, 093.B- 
0236, 094.B-0181, 095.B-0046, 096.B-0162, 097.B-0304, 
098.B-0030, 099.B-0034, 0100.B-0131, 0101.B-0156, 0102.B- 
0174, 0103.B-0165, 0104.B-0081, 0106.20ZA.001 (Omega- 
Cam), 0106.20ZA.002 (FORS2), 0108.21 YF; and under INAF 
programs 110.256C, 112.266Q; and (ii) the Liverpool Telescope, 
which is operated on the island of La Palma by Liverpool John 
Moores University in the Spanish Observatorio del Roque de 
los Muchachos of the Instituto de Astrofísica de Canarias with 
financial support from the United Kingdom Science and Tech- 
nology Facilities Council, and (iii) telescopes of the Las Cum- 
bres Observatory Global Telescope Network. 

In addition to the currently active DPAC (and ESA science) 
authors of the peer-reviewed papers accompanying Gaia DR3, 
there are large numbers of former DPAC members who made 
significant contributions to the (preparations of the) data process- 
ing. Among those are, in alphabetical order: Stephanie Accart, 
Christopher Agard, Juan José Aguado, Michaél Ajaj, Fernando 
Aldea-Montero, Alexandra Alecu, Bruno Alessi, Peter Allan, 
France Allard, Walter Allasia, Carlos Allende Prieto, Javier 
Alvarez Cid-Fuentes, Marco Antonio Alvarez, João Alves, 
Antonio Amorim, Kader Amsif, Alexandre Andrei, Antonino 
Angi, Guillem Anglada-Escudé, Erika Antiche, Sonia Anton, 
Bernardino Arcay, Clément Arnaudeau, Borja Arroyo Gal- 
ende, Vladan Arsenijevic, Tri Astraatmadja, Rajesh Kumar 
Bachchan, Adrien Bangma, Carlos Barata, Domenico Bar- 
bato, Fabio Barblan, Paul Barklem, Mickael Batailler, Duncan 
Bates, Alexandre Baudesson-Stella, Mathias Beck, Luigi Bedin, 
Dan Beilis, Antonio Bello Garcia, Vasily Belokurov, Philippe 
Bendjoya, Angel Berihuete, Hans Bernstein’, Olivier Bien- 
aymé, Lionel Bigot, Albert Bijaoui, Louis Bil, Françoise 
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Billebaud, Nadejda Blagorodnova, Thierry Bloch, Klaas de 
Boer’, Marco Bonfigli, Giuseppe Bono, Simon Borgniet, Raul 
Borrachero-Sanchez, François Bouchy, Steve Boudreault, Geral- 
dine Bourda, Guy Boutonnet, Lorenzo Bramante, Pascal Branet, 
Maarten Breddels, Scott Brown, Pierre-Marie Brunet, Thomas 
Brüsemeister, Peter Bunclark!, Roberto Buonanno, Alexandru 
Burlacu, Robert Butorafuchs, Joan Cambras, Heather Campbell, 
Sylvain Cannizzo, Christophe Carret, Manuel Carrillo, César 
Carrión, Pau Castro Sampol, Francisco Javier Casquero, Lau- 
rence Chaoul, Jonathan Charnas, Fabien Chéreau, Vincenzo 
Chiaramida Mathurin Chritin, Maria-Rosa Cioni, Uma Cladel- 
las Sanjuan, Marcial Clotet, Gabriele Cocozza, Ross Collins, 
Gabriele Comoretto, Gabriele Contursi, Leonardo Corcione, 
Gráinne Costigan, Francoise Crifo, Alessandro Crisafi, Nick 
Cross, Jan Сиурегѕ', Jean-Charles Damery, Anastasios Daper- 
golas, Eric Darmigny, Pedro David, Jonas Debosscher, Peter 
De Cat, Domitilla De Martino, Rafael De Souza, Enrique Del 
Pozo, Héctor Delgado, David Delhoume, Céline Delle Luche, 
Markus Demleitner, Léo Denglos, Sékou Diakite, Paola Di Mat- 
teo, Carla Domingues, Sandra Dos Anjos, Laurent Douchy, 
Petros Drazinos, Pierre Dubath, Javier Durán, Yifat Dzigan, 
Bengt Edvardsson, Deepak Eappachen, Sebastian Els, Arjen van 
Elteren, Kjell Eriksson, Pilar Esquej, Carolina von Essen, Wyn 
Evans, Guillaume Eynard Bontemps, Antonio Falcáo, Martí 
Farràs Casas, Jacopo Federici, Luciana Federici, Fernando de 
Felice, Agnés Fienga, Krzysztof Findeisen, Christian Fischer, 
Florin Fodor, Yori Fournier, Frédéric Franke, Benoit Frezouls, 
Aidan Fries, Jan Fuchs, Flavio Fusi Pecci, Diego Fustes, Dun- 
can Fyfe, Eva Gallardo, Silvia Galleti, Fernando García, Alberto 
García Gutiérrez, María García-Reinaldos, Daniele Gardiol, 
Nora Garralda Torres, Emilien Gaudin, Alvin Gavel, Marwan 
Gebran, Yoann Gérard, Nathalie Gerbier, Joris Gerssen, Miguel 
Gomes, Roy Gomel, Anita Gómez, Ana González-Marcos, Juan 
González-Nüfiez, Juan José González-Vidal, Eva Grebel, Michel 
Grenon, Bjórn Grieger, Eric Grux, Alain Gueguen, Pierre Guill- 
out, Julie Guiraud, Andrés Gürpide, Leanne Guy, Jean-Louis 
Halbwachs, Marcus Hauser, Aurelien Hees, Kevin Henares, 
Julien Heu, Albert Heyrovsky, Thomas Hilger, Nathan Himpens, 
Natalia Hladczuk, Wilfried Hofmann, Erik Høg, David Hogg, 
Andrew Holland, Greg Holland, Gordon Hopkinson', Claude 
Huc, Pablo Huijse, Jason Hunt, Brigitte Huynh, Arkadiusz 
Hypki, Giacinto Iannicola, Sergio Ibarmia, Vilma Icardi, Laura 
Inno, Mike Irwin, Yago Isasi Parache, Javier Izquierdo, Maja 
Jabłońska, Thierry Jacq, Asif Jan, Anne-Marie Janotto, Kevin 
Jardine, Gérard Jasniewicz, Anne Jean-Antoine Piccolo, Lau- 
rent Jean-Rigaud, Isabelle Jégouzo-Giroux, Christian Jezequel, 
Frangois Jocteur-Monrozier, Paula Jofré, Anthony Jonckheere, 
Peter Jonker, Áron Juhász, Francesc Julbe, Antonios Karam- 
pelas, Lea Karbevska, Ralf Keil, Adam Kewley, Dae-Won Kim, 
Peter Klagyivik, Jochen Klar, Jonas Klüter, Jens Knude, Angela 
Kochoska, Oleg Kochukhov, Katrien Kolenberg, Indrek Kolka, 
Pavel Koubsky, Janez Kos, Irina Kovalenko, Daniel Krefl, Maria 
Kudryashova, Ilya Kull, Alex Kutka, Frédéric Lacoste-Seris, 
Sylvain Lafosse, Valéry Lainey, Pascal Laizeau, Yannick Lasne, 
Antoni Latorre, Felix Lauwaert, Claudia Lavalley, Jean-Baptiste 
Lavigne, David Le Bouquin, Jean-Frangois Le Campion, 
Isabelle Lecoeur-Taibi, Yann Le Fustec, Vassili Lemaitre, Hel- 
mut Lenhardt, Christophe Le Poncin-Lafitte, Frédéric Leroux, 
Thierry Levoir, Hans Lindstróm, Tim Lister, Chao Liu, Mauro 
López Del Fresno, Davide Loreggia, Denise Lorenz, Cristina 
Luengo, Ian MacDonald, Marc Madaule, Pau Madrero Pardo, 
Tiago Magalhães Fernandes, Arrate Magdaleno Romeo, Kir- 
ill Makan, Valeri Makarov, Jean-Christophe Malapert, Sandrine 
Managau, Hervé Manche, Carmelo Manetta, Gregory Mantelet, 


José Marcos, Miguel Marcos Santos, Federico Marocco, Gabor 
Marschalko, Mathieu Marseille, Christophe Martayan, Óscar 
Martínez-Rubi, Michele Martino, Paul Marty, Nicolas Mary, 
Davide Massari, Benjamin Massart, Gal Matijevič, Mohamed 
Meharga, Emmanuel Mercier, Maria Messineo, Frédéric Mey- 
nadier, Daniel Michalik, Anthony Michon, Shan Mignot, Hadi 
Minbashian, Bruno Miranda, László Molnár, Marco Moli- 
naro, Giacomo Monari, Marc Moniez, Ángel Montero, Alain 
Montmory, Roger Mor, Thierry Morel, Stephan Morgenthaler, 
Angelo Mulone, Ulisse Munari, Daniel Muñoz, Cillian Mur- 
phy, Jéróme Narbonne, Gijs Nelemans, Anne-Thérése Nguyen, 
Luciano Nicastro, Sara Nieto, Thomas Nordlander, Alexandre 
Nouvel, Louis Noval, Markus Nullmeier, Derek O'Callaghan, 
Francisco Ocaña, Pierre Ocvirk, Alex Ogden, Joaquín Ordieres- 
Meré, Diego Ordonez, Giuseppe Orrü, Patricio Ortiz, José 
Osinde, Jose Osorio, Dagmara Oszkiewicz, Alex Ouzounis, 
Hugo Palacin, Max Palmer, Aviad Panahi, Chantal Panem, Vin- 
cent Papy, Peregrine Park, Ester Pasquato, Xavier Passot, Stefan 
Payne-Wardenaar, Louis Pegoraro, Roselyne Pedrosa, Christian 
Peltzer, Hanna Pentikäinen, Xavier Peñalosa Esteller, Jordi Per- 
alta, Rubén Pérez, Jean-Marc Petit, Fabien Péturaud, Bernard 
Pichon, Tuomo Pieniluoma, Anna Marina Piersimoni, François- 
Xavier Pineau, Enrico Pigozzi, Federic Pireddu, Bertrand Plez, 
Joel Poels’, Aurelian Polidoro, Eric Poujoulet, Arnaud Poulain, 
Guylaine Prat, Thibaut Prod’ homme, Andrej Prša, Elena Racero, 
Adrien Raffy, Silvia Ragaini, Serena Rago, Nicolas Rambaux, 
Piero Ranalli, Gregor Rauw, Andrew Read, José Rebordao, 
Philippe Redon, Rita Ribeiro, Ariadna Ribes Metidieri, Pas- 
cal Richard, Phil Richards, Carlos Ríos Díaz, Daniel Risquez, 
Adrien Rivard, Clement Robin, Brigitte Rocca-Volmerange, 
Maroussia Roelens, Hervé Rogues, Laurent Rohrbasser, Nico- 
las de Roll, Julia Roquette, Siv Rosén, Frederic Royer, Ste- 
fano Rubele, Laura Ruiz Dern, Idoia Ruiz-Fuertes, Federico 
Russo, Jan Rybizki, Albert Sáez Núñez, Jesús Salgado, Euge- 
nio Salguero, Nik Samaras, Paula Sánchez Gayet, Víctor 
Sánchez Giménez, Toni Santana, Helder Savietto, Maud Segol, 
Juan Carlos Segovia, Damien Segransan, Léa Sellahannadi, 
Didier Semeux, I-Chun Shih, Hassan Siddiqui, Lauri Siltala, 
André Silva, Helder Silva, Arturo Silvelo, Dimitris Sinachopou- 
los, Christos Siopis, Riccardo Smareglia, Kester Smith, Michael 
Soffel, Sergio Soria Nieto, Danuta Sosnowska, Alessandro 
Spagna, Maxime Spano, Lorenzo Spina, Ulrike Stampa, Craig 
Stephenson, Hristo Stoev, Vytautas Straizys, Frank Suess, Maria 
Siiveges, Elza Szegedi-Elek, Francis Tache, Jeff Tambouez, 
Guy Tauran, Dirk Terrell, David Terrett, Pierre Teyssandier, 
Stephan Theil, William Thuillot, Carola Tiede, Brandon Tin- 
gley, Krešimir Tisanié, Anastasia Titarenko, Jordi Torra!, Scott 
Trager, Licia Troisi, Paraskevi Tsalmantza, David Tur, Stefano 
Uzzi, Mattia Vaccari, Frédéric Vachier, Emmanouil Vachlas, 
Marc Vaillant, Gaetano Valentini, Pau Vallés, Veronique Valette, 
Emmanuel van Dillen, Walter Van Hamme, Eric Van Hemelryck, 
Wouter van Reeven, Mihaly Varadi, Marco Vaschetto, Jovan Vel- 
janoski, Lionel Veltz, Sjoert van Velzen, Teresa Via, Yves Viala, 
Jenni Virtanen, Antonio Volpicelli, Holger Voss, Viktor Votruba, 
Stelios Voutsinas, Jean-Marie Wallut, Gavin Walmsley, Olivier 
Wertz, Thomas Wevers, Rainer Wichmann, Mark Wilkinson, 
Abdullah Yoldas, Patrick Y vard, Petar Zečević, Tim de Zeeuw, 
Maruska Zerjal, Houri Ziaeepour, Claude Zurbach, and Sven 
Zschocke. 

In addition to the DPAC consortium, past and present, there 
are numerous people, mostly in ESA and in industry, who have 
made or continue to make essential contributions to Gaia, for 
instance those employed in science and mission operations or in 
the design, manufacturing, integration, and testing of the space- 
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craft and its modules, subsystems, and units. Many of those will 
remain unnamed yet spent countless hours, occasionally during 
nights, weekends, and public holidays, in cold offices and dark 
clean rooms. At the risk of being incomplete, we specifically 
acknowledge, in alphabetical order, from Airbus DS (Toulouse): 
Alexandre Affre, Marie-Thérése Aimé, Audrey Albert, Aurélien 
Albert-Aguilar, Jeanine Alloun-Etcheberry, Hania Arsalane, 
Arnaud Aurousseau, Denis Bassi, Franck Bayle, Bernard Bayol, 
Pierre-Luc Bazin, Emmanuelle Benninger, Philippe Bertrand, 
Jean-Bernard Biau, François Binter, Cédric Blanc, Eric Blonde, 
Patrick Bonzom, Bernard Bories, Jean-Jacques Bouisset, Joél 
Boyadjian, Isabelle Brault, Corinne Buge, Bertrand Calvel, Jean- 
Michel Camus, France Canton, Lionel Carminati, Michel Car- 
rie, Didier Castel, Philippe Charvet, François Chassat, Fab- 
rice Cherouat, Ludovic Chirouze, Michel Choquet, Claude 
Coatantiec, Emmanuel Collados, Philippe Corberand, Christelle 
Dauga, Robert Davancens, Catherine Deblock, Eric Decourbey, 
Charles Dekhtiar, Michel Delannoy, Michel Delgado, Damien 
Delmas, Emilie Demange, Victor Depeyre, Isabelle Desenclos, 
Christian Dio, Kevin Downes, Marie-Ange Duro, Eric Ecale, 
Omar Emam, Elizabeth Estrada, Coralie Falgayrac, Benjamin 
Farcot, Claude Faubert, Frédéric Faye, Sébastien Finana, Gré- 
gory Flandin, Loic Floury, Gilles Fongy, Michel Fruit, Flo- 
rence Fusero, Christophe Gabilan, Jérémie Gaboriaud, Cyril 
Gallard, Damien Galy, Benjamin Gandon, Patrick Gareth, Eric 
Gelis, André Gellon, Laurent Georges, Philippe-Marie Gomez, 
José Goncalves, Frédéric Guedes, Vincent Guillemier, Thomas 
Guilpain, Stéphane Halbout, Marie Hanne, Grégory Hazera, 
Daniel Herbin, Tommy Hercher, Claude Hoarau le Papillon, 
Matthias Holz, Philippe Humbert, Sophie Jallade, Grégory Jon- 
niaux, Frédéric Juillard, Philippe Jung, Charles Koeck, Marc 
Labaysse, Réné Laborde, Anouk Laborie, Jéróme Lacoste- 
Barutel, Baptiste Laynet, Virginie Le Gall, Julien L'Hermitte, 
Marc Le Roy, Christian Lebranchu, Didier Lebreton, Patrick 
Lelong, Jean-Luc Leon, Stephan Leppke, Franck Levallois, 
Philippe Lingot, Laurant Lobo, Didier Loche, Céline Lopez, 
Jean-Michel Loupias, Carlos Luque, Sébastien Maes, Bruno 
Mamdy, Denis Marchais, Alexandre Marson, Benjamin Massart, 
Rémi Mauriac, Philippe Mayo, Caroline Meisse, Hervé Mer- 
cereau, Olivier Michel, Florent Minaire, Xavier Moisson, David 
Monteiro, Denis Montperrus, Boris Niel, Cédric Papot, Jean- 
Frangois Pasquier, Gareth Patrick, Pascal Paulet, Martin Pec- 
cia, Sylvie Peden, Sonia Penalva, Michel Pendaries, Philippe 
Peres, Grégory Personne, Dominique Pierot, Jean-Marc Pillot, 
Lydie Pinel, Fabien Piquemal, Vincent Poinsignon, Maxime 
Pomelec, André Porras, Pierre Pouny, Severin Provost, Sébastien 
Ramos, Fabienne Raux, Audrey Rehby, Florian Reuscher, 
Xavier Richard, Nicolas Riguet, Mickael Roche, Gilles Rougier, 
Bruno Rouzier, Stephane Roy, Jean-Paul Ruffie, Frédéric Safa, 
Heloise Scheer, Claudie Serris, André Sobeczko, Jean-François 
Soucaille, Romain Suze, Philippe Tatry, Théo Thomas, Pierre 
Thoral, Dominique Torcheux, Vincent Tortel, Damien Tourbez, 
Stephane Touzeau, Didier Trantoul, Cyril Vétel, Jean-Axel 
Vatinel, Jean-Paul Vormus, and Marc Zanoni; from Airbus DS 
(Friedrichshafen): Jan Beck, Frank Blender, Volker Hashagen, 
Armin Hauser, Bastian Hell, Cosmas Heller, Matthias Holz, 
Heinz-Dieter Junginger, Klaus-Peter Koeble, Karin Pietroboni, 
Ulrich. Rauscher, Rebekka Reichle, Florian Reuscher, Ari- 
ane Stephan, Christian Stierle, Riccardo Vascotto, Christian 
Hehr, Markus Schelkle, Rudi Kerner, Udo Schuhmacher, Peter 
Moeller, Rene Stritter, Jürgen Frank, Wolfram Beckert, Eve- 
lyn Walser, Steffen Roetzer, Fritz Vogel, and Friedbert Zilly; 
from Airbus DS (Stevenage): Mohammed Ali, Bill Bental, David 
Bibby, Leisha Carratt, Veronica Carroll, Clive Catley, Patrick 
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Chapman, Christoper Chetwood, Alison Colegrove, Tom Cole- 
grove, Andrew Davies, Denis Di Filippantonio, Andy Dyne, 
Alex Elliot, Omar Emam, Colin Farmer, Steve Farrington, Nick 
Francis, Albert Gilchrist, Brian Grainger, Yann Le Hiress, Vicky 
Hodges, Jonathan Holroyd, Haroon Hussain, Roger Jarvis, 
Lewis Jenner, Steve King, Chris Lloyd, Neil Kimbrey, Alessan- 
dro Martis, Bal Matharu, Karen May, Florent Minaire, Kather- 
ine Mills, James Myatt, Chris Nicholas, Paul Norridge, David 
Perkins, Michael Pieri, Matthew Pigg, Angelo Povoleri, Robert 
Purvinskis, Phil Robson, Julien Saliege, Satti Sangha, Paramijt 
Singh, John Standing, Dongyao Tan, Keith Thomas, Rosalind 
Warren, Andy Whitehouse, Robert Wilson, Hazel Wood, Steven 
Danes, Scott Englefield, Juan Flores-Watson, Chris Lord, Allan 
Parry, Juliet Morris, Nick Gregory, and Ian Mansell. 

From ESA, in alphabetical order: Ricard Abello, Asier 
Abreu, Ivan Aksenov, Matthew Allen, Salim Ansari, Philippe 
Armbruster, Mari-Liis Aru, Alessandro Atzei, Liesse Ayache, 
Samy Azaz, Nana Bach, Jean-Pierre Balley, Paul Balm, Manuela 
Baroni, Rainer Bauske, Thomas Beck, Gabriele Bellei, Car- 
los Bielsa, Gerhard Billig, Carmen Blasco, Andreas Boosz, 
Bruno Bras, Julia Braun, Thierry Bru, Frank Budnik, Joe 
Bush, Marco Butkovic, Jacques Candeé, David Cano, Carlos 
Casas, Francesco Castellini, David Chapmann, Nebil Cinar, 
Mark Clements, Giovanni Colangelo, Peter Collins, Ana Col- 
orado McEvoy, Gabriele Comoretto, Vincente Companys, Fed- 
erico Cordero, Yannis Croizat, Sylvain Damiani, Fabienne 
Delhaise, Gianpiero Di Girolamo, Yannis Diamantidis, John 
Dodsworth, Ernesto Dólling, Jane Douglas, Jean Doutreleau, 
Dominic Doyle, Mark Drapes, Frank Dreger, Peter Droll, 
Gerhard Drolshagen, Michal Ďurovič, Bret Durrett, Christina 
Eilers, Yannick Enginger, Alessandro Ercolani, Matthias 
Erdmann, Orcun Ergincan, Robert Ernst, Daniel Escolar, Maria 
Espina, Hugh Evans, Fabio Favata, Stefano Ferreri, Daniel 
Firre, Michael Flegel, Melanie Flentge, Alan Flowers, Steve 
Foley, Julia Fortuno-Benavent, Jens Freihöfer, Rob Furnell, 
Julio Gallegos, Maria De La Cruz Garcia Gonzalez, Philippe 
Garé, Wahida Gasti, José Gavira, Frank Geerling, Franck Ger- 
mes, Gottlob Gienger, Bénédicte Girouart, Bernard Godard, 
Nick Godfrey, César Gómez Hernández, Roy Gouka, Cosimo 
Greco, Robert Guilanya, Kester Habermann, Manfred Had- 
wiger, Ian Harrison, Angela Head, Martin Hechler, Javier Her- 
nando Bravo, Kjeld Hjortnaes, John Hoar, Jacolien Hoek, 
Frank Hoffmann, Justin Howard, Fredrik Hülphers, Arjan 
Hulsbosch, Christopher Hunter, Premysl Janik, José Jiménez, 
Emmanuel Joliet, Helma van de Kamp-Glasbergen, Simon 
Kellett, Andrea Kerruish, Kevin Kewin, Oliver Kiddle, Sabine 
Kielbassa, Volker Kirschner, Kees van 't Klooster, Ralf Koh- 
ley, Jan Kolmas, Oliver El Korashy, Arek Kowalczyk, Holger 
Krag, Benoit Lainé, Markus Landgraf, Sven Landstróm, Math- 
ias Lauer, Robert Launer, Laurence Tu-Mai Levan, Mark ter Lin- 
den, Santiago Llorente, Tim Lock’, Alejandro Lopez-Lozano, 
Guillermo Lorenzo, Tiago Loureiro, James Madison, Juan 
Manuel Garcia, Federico di Marco, Jonas Marie, Filip Marinic, 
Pier Mario Besso, Arturo Martín Polegre, Ander Martínez, 
Monica Martínez Fernández, Marco Massaro, Paolo de Meo, 
Ana Mestre, Claudio Mevi, Luca Michienzi, David Milligan, 
Ali Mohammadzadeh, David Monteiro, Richard Morgan-Owen, 
Trevor Morley, Prisca Mühlmann, Jana Mulacova, Michael 
Müller, Pablo Mufioz, Petteri Nieminen, Alfred Nillies, Wilfried 
Nzoubou, Alistair O'Connell, Karen O'Flaherty, Alfonso Olias 
Sanz, William O'Mullane, José Osinde, Oscar Pace, Mohini 
Parameswaran, Ramon Pardo, Taniya Parikh, Paul Parsons, 
Panos Partheniou, Torgeir Paulsen, Dario Pellegrinetti, José- 
Louis Pellon-Bailon, Joe Pereira, Michael Perryman, Christian 
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Philippe, Alex Popescu, Frédéric Raison, Riccardo Rampini, 
Florian Renk, Alfonso Rivero, Andrew Robson, Gerd Róssling, 
Martina Rossmann, Markus Riickert, Andreas Rudolph, Frédéric 
Safa, Johannes Sahlmann, Eugenio Salguero, Jamie Salt, Gio- 
vanni Santin, Fabio de Santis, Rui Santos, Giuseppe Sarri, 
Stefano Scaglioni, Melanie Schabe, Dominic Schafer, Micha 
Schmidt, Rudolf Schmidt, Ared Schnorhk, Klaus-Jiirgen Schulz, 
Jean Schütz, Julia Schwartz, Andreas Scior, Jórg Seifert, 
Christopher Semprimoschnig!, Ed Serpell, Iñaki Serraller Viz- 
caino, Gunther Sessler, Felicity Sheasby, Alex Short, Hassan 
Siddiqui, Heike Sillack, Swamy Siram, Chloe Sivac, Christo- 
pher Smith, Claudio Sollazzo, Steven Straw, Daniel Tapiador, 
Pilar de Teodoro, Mark Thompson, Giulio Tonelloto, Felice 
Torelli, Raffaele Tosellini, Cecil Tranquille, Irren Tsu-Silva, 
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